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ABSTRACT 

We describe a new method to achieve point spread function (PSF) subtractions for high- 
contrast imaging using Principal Component Analysis (PCA) that is applicable to both point 
sources or extended objects (disks). Assuming a library of reference PSFs, a Karhunen-Loeve 
transform of theses references is used to create an orthogonal basis of eigenimages, on which the 
science target is projected. For detection this approach provides comparable suppression to the 
Locally Optimized Combination of Images (LOCI) algorithm, albeit with increased robustness 
to the algorithm parameters and speed enhancement. For characterization of detected sources 
the method enables forward modeling of astrophysical sources. This alleviates the biases in the 
astrometry and photometry of discovered faint sources, which are usually associated with LOCI- 
based PSF subtractions schemes. We illustrate the algorithm performance using archival Hubble 
Space Telescope (HST) images, but the approach may also be considered for ground-based data 
acquired with Angular Differential Imaging (ADI) or integral-field spectrographs (IFS). 

Subject headings: methods: data analysis — techniques: image processing 



Introduction 



Imaging faint exoplanets and circumstellar disks is a difficult task limited by the small angular separation 
and h igh contrast between the faint astrophysical signal and the residual starlight (jOppenheimer fc Hinklev 



20091 ) . The first direct images of planets or other interesting localized structures orbiting nearby stars were 
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Lagrange et al. 


2009, 


2010; 


Marois et al. 



pand in the near future with a new generation ground-based instr uments being implemented ([Macintosh et al 



2008 ; Beuzit et al.ll2008t Hinklev et alj|2011 ; Hodapp et alJbooil) . Even using advan ced starlight suppression 



systems (e.g. coronag r aphs) such as the one d esigned for the Gemini Planet Imager (jSivaramakrishnan et al 



2001; Soummer 2005: Soummer et al 



2011b), image post-processing remains necessary to achieve the full 



high-contrast performance ( Marois et al. 2008bh . Indeed, most recent direct imaging results were enable d by 
key advances in observational procedures, e.g. Angular Differential Imaging (ADD (IMarois et al.1 120061). s i- 
multaneous imaging of contiguous wavelengths ( Marois et al. 2000 ; Sparks &; Ford 20021 : Puevo et al. 2012 ). 
in conjunction with sophis ticated images analysi s routines based on the Locally Optimized Combination of 
Images (LOCI) algorithm (jLafreniere et al.ll2007j) . Although LOCI is extremely powerful at detecting point 
sources, the very aggressive PS F subtraction introd uces biases on the photometry and astrometry, which 
have to be calibrated carefully (jMarois et al.l 2010a). Biases can b e reduced significantly using statistical 
exploration of the algorithm parameter space (jSoummer et al.ll2011al ) , a computationally e xpensive solution, 
or by modifyin g the underlying co st function by masking some pixels (jMarois et al.ll2010al ) or using damping 
penalty terms (jPuevo et al.l 120121 ). LOCI tends to subtract extended structures such as debris disks, and 
several studies have focused on mitigation stra tegies, taking advantage of the geometry in edge-on disks 



or pe rforming an optimization over large zones ([Fitzgerald et al.l 120071 : iThalmann et alJl2010t iBuenzli et al 
2010|) . 



In this paper we introduce a new reduction algorithm based on projections on a truncated Karhunen 
Loeve transformation of the reference PSF library. We first introduce the general formalism associated with 
PSF subtraction and describe the algorithm. We then use HST NICMOS coronagraphic data to illustrate 
its performances and highlight its advantages for the characterization of point sources and the discovery of 
faint extended structures around nearby stars. 



2. The KL Image Projection (KLIP) algorithm 

2.1. Optimal PSF subtraction 

An astronomical Point Spread Function (PSF) can be seen as a A^-dimensional stochastic process indexed 
over a scalar quantity ip = ijj{t,X,9,(p,T eX p,mstar), which is representative of the state of the instrument- 
telescope system (for instance: time, wavelength, field rotation angle, wavefront error, exposure time, star 
magnitude) at the time of the exposure. With N pixels in the image, the PSF intensity can be written as: 

7^(n) =I(p[ri\,q[n],il>), (1) 

where we have rearranged the detector coordinates in one dimension using the pixel mapping functions p 
and q and the pixel index n (n G [l,n]). 

We assume an observation sequence that yields a target image: 

T(n) = I fo (n) + eA(n), (2) 



which might (e = 
images: 



1) or might not (e = 0) contain faint astronomical signal A(n), and a set of reference 
{R k (n) = Iip k {n)} k =i...K- (3) 
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tpo (n) and ipk(n>) represent the realizations of the state of the telescope-instrument system respectively for 
the target image and the set of references. 

Using prior information about the observation strategy, these references have been selected so that they 
do not contain any astronomical signal in the search region S (dimension N$ pixels) of the PSF. Moreover the 
observations have been carried out at neighboring states of the telescope/instrument system: | tpi — 4>j \ < 5 
for all i,j € [0, K], where S is a small quantity that can be adjusted in practice using the correlation between 
images as a proxy. 

The goal of PSF subtraction is to reconstruct (n) using the set of reference images and to subtract 
it from the target image since eA{n) = T — I^ (n). Unfortunately I,p(n) is a continuous random process of 
ijj and thus an infinite amount of references would be needed to reconstruct (n) exactly. In practice we 
can only seek to determine the best estimate (n) of (n) using the limited information contained in the 
K references. 

Assuming that each reference is unique, the ensemble {Itf, k }k=i...K spans Ik, a K dimensional sub-space 
of S. Ns , which has an infinity of orthonormal bases {Zk}k=i...K- This problem of estimating I^ (n) can be 
formulated as: 



What is the basis set {Zk(n)}k=x...K 'most likely to minimize the distance between a random realization 
of lib (n) (ip = ipo in our case) andlx? More formally: 



mm 

{Z h } 




K V 

^ < I^,Zk >s Z k (n) 

k=l / 



(4) 



where E^[-] stands for the expected value over the telescope realizations and < •,• >,s denotes the inner 
product over S. This is a well known problem in signal processing and, assuming that all the images are 
of zero mean over S, the optimal basis set {Z^ L }k=i...K is given by the Karhunen-Loeve transform of the 



{Itp k }k=i...K (|Karhunenlll947[ lLoevdll948l) . Moreover if one chooses to only describe over a Kku p < K 



dimensional sub-space, then the optimal basis set is given by the truncated Karhunen-Loeve transform of 
the references, {Z£ L } k =i...K kUp 



2.2. The algorithm 

The KLIP algorithm consists of the following steps: 

1. partition the target T(n) and references Rk{n) images in an ensemble of search areas S, and subtract 
their average values so that they have zero mean over each S. 

2. compute the Karhunen-Loeve transform of the set of reference PSFs {Rk{n)}k=i...K- 

Z?\n) = ^=E c k (iP P )R P (n\ (5) 
V A k p=1 

where the vectors Ck = [ck{ipi)---c-k{' l pK )] are the eigenvectors of the K-dimensional covariance matrix 
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Err of the Rk(n) references over S, and {A}fc = i...jc are its eigenvalues. 



EjirM = J^Rii^R^n) (6) 

n=l 

^Kfl • C k = AfeCfe, A! > A 2 > .... > A N (7) 



3. choose a number of modes K^i ip to keep in the estimate J^, (n). 

4. compute the best estimate I^, (n) of the actual PSF (n) from the projection of the science target 



T(n) on the {Zjf L (ri)}k=x,,,K k H Karhunen-Loeve basis: 



Uo(n)= E <T,Z^> s Z^ L (n). (8) 
fc=i 



5. calculate the final image F(n) = T(n) — j^, (n), 



F(n) = ( J^(n) -J2< Uo,Zk L >s Z« L (n) 
fc=i 



6 4(n) - E < A L > 5 » . (9) 



fc=i 



2.3. Throughput and Signal to Noise Ratio 

Equation|9]provides fundamental information pertaining to the optimal choice of Kku p for a given search 
zone iS. The variance in the search zone S after the PSF subtraction is given by: 

Kuip / „ T ry KL _ x 2 



fe=l ^ 



A t 



(11) 



which is a decreasing function of Kku p . In Figure [T] we show the impact of Kku p on the residual variance 
integrated over entire KLIP-reduced images (i.e. S is the full image), using HST-NICMOS coronagraphic 
data preselected to have no astrophysical signal. Figure [T] shows that the empirical residual variance after 
PSF subtraction follows the theoretical estimate in Equation [TT1 

Similarly, equation |H] provides a simple expression for the signal: 

A{n)- <A,Z^ L > s Z^ L (n) j ,12) 

k=l 

For the dominant modes (i.e. k <C K) < A, Z¥ L >g ~ since the morphology of the astronomical signal 
cannot be reproduced by a PSF realization. Indeed the PSF of a faint planet is localized over a few pixels 
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and thus almost orthogonal to the main modes of the telescope's PSF realizations. Naturally this argument 
depends on the spatial morphology of the signal: in Figure [5] we show that the algorithm throughput 
decreases when Kku p increases for several types of faint astrophysical object. Given priors on the shape of 
the astrophysical signal (i.e. planet position and flux, or a disk model), it is possible to estimate directly the 
Signal to Noise Ratio (SNR) as a function Kku p when S is small enough so that equation [TT] represents the 
local noise in the vicinity of the source. It is therefore possible to choose K^up a priori in order to optimize 
the SNR for a given putative source in the image. 



i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — r 




50 100 150 200 250 300 

KLIP truncation 

Fig. 1. — Normalized total residual variance over the entire image after PSF subtraction. The black curve 
corresponds to the theoretical expression cr| / 1 1 /,/,„ 1 1 2 based on the corollary of the Karhunen-Loeve theorem 
(See Equation I lip . The color curves correspond to the empirical total residual variance after subtraction for 
several target stars from the HST NICMOS archive without any astrophysical signal around them. In the 
case of a target with an astrophysical signal, the variance saturates to a minimum level and does not follow 
the theoretical expression. The theoretical expression from Equation [TT] can therefore be used to estimate 
semi-analytically the residual noise after subtraction. If KLIP is applied in small regions (e.g. using zones as 
with LOCI) this expression can be used to estimate the noise locally around a putative source in the target 
image. 



3. Discussion 

3.1. Forward Modeling with KLIP 

In practice the optimal search area S might be too large to use Equation [TTJ and Equation [12] to pre-select 
Kkup- Instead Kku p can be tuned by inspecting the SNR of putative sources in reduced images. However 
the second term of Equation [5] carries all the information necessary to characterize faint astrophysical signal 
once it has been detected. Indeed the impact of the algorithm on the signal only depends on the projection of 
A(n) on the {Zjf L }k=i..K klip and is not a function of I^ . As a consequence KLIP enables forward modeling 
of detected planets and disks by fitting directly an astrophysical model A^(n) to the final reduced images 
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KLIP truncation 



Fig. 2. — KLIP throughput based on the expression in Equation [9] from the projection of the astrophysical 
signal A onto the KL basis < A, Zjf L >s- Since a planet's signal is mostly orthogonal to the main KL modes 
the throughput is high, as shown for a planet located at ~ 1.5 arcsec from the star (top curve). For more 
extended sources the throughput is reduced, as illustrated for a simple geometrical model corresponding to 
the HD 181327 disk (middle curve). The bottom curve corresponds to a hypothetical ring-like disk about 
three times wider than HD 181327. Note that here iS corresponds to an entire HST NICMOS image, the 
throughput will be smaller when using smaller partitions of the image. However because it is calculated 
analytically for a source model independent from its flux, it can be calibrated easily. 
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F(ri), where £ represents the model parameters (e.g. astrometry and photometry for a planet, or geometry, 
grain distribution for a disk). £ can be obtained by minimizing the following cost function: 

Kklip 



min < 



£ F(n) A e (n) + £ < A c , Zf * > 5 Z** 

n=l \ fc=l 




(13) 



A similar approach was proposed bv lMarois et al.l (|2010a[) using LOCI, but it presents degeneracies because 
the LOCI coefficients are not guaranteed to be independent of or A. Such degeneracies do not arise in 
the case of KLIP because the effect of the algorithm on the astrophysical source is completely determined 
from Equation [TBI and because the basis i Z^ L (n)}k = ^ ...Kk u* ^ s independent of the science target. Note that 



Eg uat ion IT31 assumes a noiseless signal. Optimal observers ( Kasdin fc Braems]|2006 ; Caucci et al 1 l2007il2009h 
can mitigate the impact of noise in the final estimate of £, and their implementation is greatly facilitated by 
the fact that Equation[T3lonlv consists of propagating A^(n) through a linear filter. Discussion regarding the 
practical implementation and performances of such observers are however beyond the scope of this study. 

The left panel of Figure [3] shows an archival HST NICMOS image of HR 8799 reduced with KLIP using 
annular search zon es. While the th r ee plan ets 6,c, and d are visible they are not detected unambiguously: the 
result reported by ISoummer et al.l (|2011al) used a combination of thousands of LOCI reductions, and took 
advantage of two HST roll angles to minimize the false alarm probability. One could envision a similar pa- 
rameter search with KLIP. The computational cost of such an exercise would howeve r be reduced since KLIP 
does not require the introduction of "optimization zones" and "subtraction zones" (jLafreniere et al J 12007). 
which reduces the dimensionali t y of th e algorithm's parameter exploration. More importantly the astrome- 
try reported in ISoummer et al.l (|2011ah was the result of a lengthy characterization of the biases associated 
with LOCI, involving thousands of reductions with synthetic companions. On the right panel of Figure [3] 
we illustrated how one can determine the astrometric location and the photometry of HR8799bcd with for- 
ward modeling in a single KLIP-reduced image. In this example we performed a simple \ 2 minimization of 
the position and flux of each planets, but this framework enables easy implementation of more appropriate 
modeling using Markov Chains Monte Carlo (MCMC). KLIP thus provides astrophysical estimates which 
do not contain any bias arising from self-subtraction in the reduction algorithm. The uncertainties on the 
astrophysical parameters are solely driven by the local noise in the reduced image. This can be addressed 
by optimizing the KLIP parameters to increase local SNR, using several exposures and taking advantage of 
observers that are optimal with respect to the residual noise. 



3.2. Similarities between LOCI and KLIP 



PCA ha s been extensively used in a stronomy to reduce the dimensionality of large datasets, e.g. 
Cowan et al. ( 20091 ): Connolly et al. (1995). Principal components can also be used directly with LOCI 
or any of its variants to reduce the dimensionality of the reference set, using KL eigen- images instead of 
reference PSFs. 

The co st function associated with KLIP (Equation [4| is very similar to the cost function associated 
with LOCI (jLafreniere et alj 120071 ). However, LOCI minimizes the distance between the science target and 
the reference PSFs, whereas KLIP minimizes, in the statistical sense, the distance between any PSFs 1^, in 
the close vicinity of references library and the i^jp-dimensional orthonormal basis-set {Zj? L (n)}k=i...K k n p - 
It is in principle possible to include the target to the reference library before decomposition into principal 
components, but in this case astrophysical sources can be completely subtracted and forward modeling cannot 



Fig. 3.— Left: Final KLIP-subtracted image from HR8799 using HST NICMOS archival data from 1998. 
The result is obtained by applying KLIP locally in large annular zones using the first 42 KL images. Right: 
result from forward modeling using Equation [T3J Note that the left image correspo nds to a single reductio n 
and is comparable in contrast to single images obtained previously with LOCI (jSoummer et al.l l2011a[) . 
Deeper images can be obtained by combining data from different orients and by exploring the algorithm 
parameter space. 



be carried out. While such solutions might yield slightly deeper contrasts they might yield astrophysical 
estimates with biases that are difficult to calibrate. In this sense LOCI, with well chosen parameters, might 
be optimal for detection while KLIP is more robust for characterization. 

Because the set of {i?fc}fc=i...jv is not an orthonormal basis, LOCI requires the inversion of the covariance 
matrix Err, which is typically ill-conditioned. This leads to noise amplification if the inversion i s not properly 



regularized. Several approaches to regularization exist, for example eigenvalue truncation (Marois et al 



2010a ) , or Tikhonov regularization ( Puevo et al.ll2012 ; Soummer et al 1 l2011ah . With KLIP the regularizat 



tion 

is implicit with the truncation of the KL basis. It is also more easily controlled because the KL set remains 
optimal for any value of Kku p ■ In addition, simple regularization schemes do not reduce the dimensionality of 
the covariance matrix, whereas KLIP reduces the number of references and therefore improves speed. Once 
the KL transform has been generated for a given search area geometry, the remaining free parameter Kku p 
can be explored with minimal computation burden (projections). Depending on the implementation of the 
algorithm and the type of data, KLIP can lead to considerable speed increases, especially when considering 
KLIP does not require to process synthetic sources for characterization of detected signal. 

As far as detection is concerned the most striking difference between LOCI and KLIP resides in disk 
imaging. Using HST Cycle 10 NICMOS archival data of the disk around HD 181327, we performed a global- 
LOCI and global-KLIP reduction over the entire image, using a carefully pre-selected reference library of 232 
PSFs. The disk edges imaged with KLI P, Fig. [H are significa ntly better constrained when compared with 



rb^s. ine disk edges imaged with rvLir, .big. |4J are significantly better constrained when compared with 
state of the art classical PSF subtraction ( Schneider et al.fl999f ). For comparison we show in Figure|4]a LOCI 



reduction of the same object without any regularization: the self subtraction degrades the rich astrophysical 
information contained in the ring's edge. It is important to note that with appropriate regularization it is 
possible to obtain with LOCI a very similar result to KLIP. However, with LOCI it is not easy to model the 
astrophysical signal since there is no guarantee that the algorithm configuration optimal for detection will 
not bias the astrophysical signal. 



Fig. 4. — left: image of HD 181327 using KLIP from a carefully selected reference library of 232 PSFs, 
using 35 KLIP coefficients. For comparison the right panel shows the image obtained with LOCI without 
regularization using the same dataset. We were able to obtain a very similar image using appropriate 
regularization with LOCI. However, there is no guarantee that the structure detected in a reduced LOCI 
image is not an artifact of the regularization. KLIP circumvents this degeneracy and enables forward 
modeling. 

4. Conclusion 

In this paper we have introduced a new algorithm for high-contrast image post-processing. KLIP is a 
cousin of existing methods (e.g LOCI) in the sense that it builds, using a quadratic cost function, a synthetic 
reference PSF as a linear combination of references from a library. KLIP involves a truncated Karhunen- 
Loeve transform of the reference library in order to produce an optimal set of orthogonal references. The 
science target is then projected onto a truncated set of these eigenimages to generate the synthetic reference. 
We used HST NICMOS data to show that KLIP provides detectability performances similar to LOCI. 
In addition to reducing the dimensionality of the data with the KL truncation, KLIP provides a linear 
framework that enables forward modeling of astrophysical sources by fitting directly an astrophysical model 
to the PSF-subtracted data, without introducing degeneracies. For example the algorithm throughput for a 
planet candidate can be calculated semi-analytically and independently of the planet signal without having 
to inject synthetic companions in a reference PSF. Taking advantage of the properties of the KL transform, 
it is possible to generate a theoretical expression of the SNR in subtracted images, provided that KLIP is 
applied is sufficiently small zones to be representative of the local noise around the planet candidate. The 
algorithm was illustrated using HST data, but it is also applicable to ground-based IFS data, and preliminary 
tests show that a local application is preferable in this case because of the lesser degree of correlation between 
PSFs. 

The authors also thank J. B. Hagan, M. Perrin, S. Lonsdale, T. Comeau and M. Russ for work and help 
on pipeline development; G. Schneider and D. Hines for discussions on disk imaging with NICMOS, and for 
providing the LAPL reference PSF library. The authors also thank C. Chen, M. Fitzgerald, B. Menard, C. 
Norman, B. Oppenheimer, A. Sivaramakrishnan, A. Szalay, G. Vasisht for discussions. Support for Program 
number HST-AR-12652.01-A was provided by NASA through a grant from the Space Telescope Science 
Institute, which is operated by the Association of Universities for Research in Astronomy, Incorporated, 



-10- 



undcr NASA contract NAS5-26555. This work was also performed in part under contract with the Jet 
Propulsion Laboratory (JPL) funded by NASA through the Sagan Fellowship Program. JPL is managed for 
NASA by the California Institute of Technology. 



REFERENCES 

Beuzit, J.-L., Feldt, M., Dohlen, K., Mouillet, D., Puget, P., Wildi, F., Abe, L., Antichi, J., Baruffolo, A., 
Baudoz, P., Boccaletti, A., Carbillet, M., Charton, J., Claudi, R., Downing, M., Fabron, O, Feautrier, 
P., Fcdrigo, E., Fusco, T., Gach, J.-L., Gratton, R., Henning, T., Hubin, N., Joos, F., Kasper, M., 
Langlois, M., Lenzen, R., Moutou, C, Pavlov, A., Petit, C, Pragt, J., Rabou, P., Rigal, F., Roelfsema, 
R., Rousset, C, Saisse, M., Schmid, H.-M., Stadler, E., Thalmann, C, Turatto, M., Udry, S., Vakili, 
F., & Waters, R. 2008, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference 
Series, Vol. 7014, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series 

Buenzli, E., Thalmann, O, Vigan, A., Boccaletti, A., Chauvin, C, Augereau, J. O, Meyer, M. R., Menard, 
F., Desidera, S., Messina, S., Henning, T., Carson, J., Montagnier, G., Beuzit, J. L., Bonavita, M., 
Eggenberger, A., Lagrange, A. M., Mesa, D., Mouillet, D., & Quanz, S. P. 2010, A&A, 524, LI 

Caucci, L., Barrett, H. H., Devaney, N., & Rodriguez, J. J. 2007, Journal of the Optical Society of America 
A, 24, 13 

Caucci, L., Barrett, H. H., & Rodriguez, J. J. 2009, Optics Express, 17, 10946 

Connolly, A. J., Szalay, A. S., Bcrshady, M. A., Kinney, A. L., & Calzetti, D. 1995, AJ, 110, 1071 

Cowan, N. B., Agol, E., Meadows, V. S., Robinson, T., Livcngood, T. A., Deming, D., Lisse, C. M., A'Hearn, 
M. F., Wellnitz, D. D., Seager, S., Charbonncau, D., & EPOXI Team. 2009, ApJ, 700, 915 

Fitzgerald, M. P., Kalas, P. G., Duchene, G., Pinte, C, & Graham, J. R. 2007, ApJ, 670, 536 

Hinkley, S., Oppenheimer, B. R., Zimmerman, N., Brenner, D., Parry, I. R., Crepp, J. R., Vasisht, G., Ligon, 
E., King, D., Soummer, R., Sivaramakrishnan, A., Beichman, C, Shao, M., Roberts, L. C, Bouchez, 
A., Dekany, R., Pueyo, L., Roberts, J. E., Lockhart, T., Zhai, C, Shelton, C, & Burruss, R. 2011, 
PASP, 123, 74 

Hodapp, K. W., Suzuki, R., Tamura, M., Abe, L., Suto, H., Kandori, R., Morino, J., Nishimura, T., Takami, 
H., Guyon, O., Jacobson, S., Stahlberger, V., Yamada, H., Shelton, R., Hashimoto, J., Tavrov, A., 
Nishikawa, J., Ukita, N., Izumiura, H., Hayashi, M., Nakajima, T., Yamada, T., & Usuda, T. 2008, 
in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 7014, Society 
of Photo-Optical Instrumentation Engineers (SPIE) Conference Series 

Janson, M., Carson, J. C, Lafreniere, D., et al. 2012, ApJ, 747, 116 

Kalas, P., Graham, J. R., Chiang, E., Fitzgerald, M. P., Clampin, M., Kite, E. S., Stapelfcldt, K., Marois, 
C, & Krist, J. 2008, Science, 322, 1345 

Kalas, P. 2011, IAU Symposium, 276, 279 

Karhunen, H. 1947, Ann. Acad. Science. Fenn, Scr. A.I. 37 

Kasdin, N. J. & Braems, I. 2006, ApJ, 646, 1260 



- 11 - 



Kraus, A. L., & Ireland, M. J. 2012, ApJ, 745, 5 

Lafreniere, D., Marois, C, Doyon, R., Nadeau, D., & Artigau, E. 2007, ApJ, 660, 770 

Lagrange, A., Bonnefoy, M., Chauvin, G., Apai, D., Ehrenreich, D., Boccaletti, A., Gratadour, D., Rouan, 
D., Mouillet, D., Lacour, S., & Kasper, M. 2010, Science, 329, 57 

Lagrange, A.-M., Gratadour, D., Chauvin, G., Fusco, T., Ehrenreich, D., Mouillet, D., Rousset, G., Rouan, 
D., Allard, F., Gendron, E., Charton, J., Mugnier, L., Rabou, P., Montri, J., & Lacombe, F. 2009, 
A&A, 493, L21 

Loeve, M. 1948, Hermann Paris, France 

Macintosh, B. A., Graham, J. R., Palmer, D. W., Doyon, R., Dunn, J., Gavel, D. T., Larkin, J., Oppenheimer, 
B., Saddlcmycr, L., Sivaramakrishnan, A., Wallace, J. K., Bauman, B., Erickson, D. A., Marois, C, 
Poyneer, L. A., & Soummer, R. 2008, in Society of Photo-Optical Instrumentation Engineers (SPIE) 
Conference Series, Vol. 7015, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference 
Series 

Marois, C, Doyon, R. ., Racine, R. ., & Nadeau, D. 2000, PASP, 112, 91 

Marois, C, Lafreniere, D., Doyon, R., Macintosh, B., & Nadeau, D. 2006, ApJ, 641, 556 

Marois, C, Macintosh, B., Barman, T., Zuckerman, B., Song, I., Patience, J., Lafreniere, D., & Doyon, R. 
2008a, Science, 322, 1348 

Marois, C, Macintosh, B., Soummer, R., Poyneer, L., & Bauman, B. 2008b, in Society of Photo-Optical 
Instrumentation Engineers (SPIE) Conference Series, Vol. 7015, Society of Photo-Optical Instrumen- 
tation Engineers (SPIE) Conference Series 

Marois, C, Macintosh, B., & Veran, J. 2010a, in Society of Photo-Optical Instrumentation Engineers (SPIE) 
Conference Series, Vol. 7736, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference 
Series 

Marois, C, Zuckerman, B., Konopacky, Q. M., Macintosh, B., & Barman, T. 2010b, Nature, 468, 1080 
Oppenheimer, B. R. & Hinkley, S. 2009, ARA&A, 47, 253 

Pueyo, L., Crepp, J. R., Vasisht, G., Brenner, D., Oppenheimer, B. R., Zimmerman, N., Hinkley, S., Parry, 
I., Beichman, C, Hillenbrand, L., Roberts, L. C, Dekany, R., Shao, M., Burruss, R., Bouchez, A., 
Roberts, J., & Soummer, R. 2012, ApJS, 199, 6 

Schneider, G., Smith, B. A., Becklin, E. E., Koerner, D. W., Meier, R., Hines, D. C, Lowrance, P. J., Terrile, 
R. J., Thompson, R. I., & Rieke, M. 1999, ApJ, 513, L127 

Sivaramakrishnan, A., Koresko, C. D., Makidon, R. B., Berkefeld, T., & Kuchner, M. J. 2001, ApJ, 552, 397 

Soummer, R. 2005, ApJ, 618, L161 

Soummer, R., Hagan, J., Pueyo, L., Thormann, A., Rajan, A., & Marois, C. 2011a, ApJ, 741, 55 
Soummer, R., Sivaramakrishnan, A., Pueyo, L., Macintosh, B., & Oppenheimer, B. R. 2011b, ApJ, 729, 144 
Sparks, W. B. & Ford, H. C. 2002, ApJ, 578, 543 



- 12 - 



Thalmann, C, Grady, C. A., Goto, M., Wisniewski, J. P., Janson, M., Henning, T., Fukagawa, M., Honda, 
M., Mulders, G. D., Min, M., Moro-Martin, A., McElwain, M. W., Hodapp, K. W., Carson, J., Abe, 
L., Brandner, W., Egner, S., Feldt, M., Fukue, T., Golota, T., Guyon, O., Hashimoto, J., Hayano, Y., 
Hayashi, M., Hayashi, S., Ishii, M., Kandori, R., Knapp, G. R., Kudo, T., Kusakabe, N., Kuzuhara, 
M., Matsuo, T., Miyama, S., Morino, J.-L, Nishimura, T., Pyo, T.-S., Scrabyn, E., Shibai, H., Suto, 
H., Suzuki, R., Takami, M., Takato, N., Terada, H., Tomono, D., Turner, E. L., Watanabe, M., 
Yamada, T., Takami, H., Usuda, T., & Tamura, M. 2010, ApJ, 718, L87 



This preprint was prepared with the AAS IATjTjX macros v5.2. 



